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Abstract 

We consider three-loop radiative-recoil corrections to hyperfine splitting in muonium due to 
insertions of one-loop polarization operator in the muon factor. The contribution produced by 
electron polarization insertions are enhanced by the large logarithm of the electron-muon mass 
ratio. We obtained all single-logarithmic and nonlogarithmic radiative-recoil corrections of order 
o?{m/M)Ep generated by the diagrams with electron and muon polarization insertions. 
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I. INTRODUCTION. MUON FACTOR CONTRIBUTION TO HFS 



Calculation of numerous corrections to hyperfme splitting (HFS) in muonium has a long 
history (see, e.g., reviews in [HI2])- By now the largest challenge to the theory is calculation 
of three-loop radiative-recoil corrections. Below we consider such corrections generated by 
two sets of diagrams with electron and muon vacuum polarization insertions in the radiative 
photon attached to the muon line in Fig. [T] and Fig. [2] The muon anomalous magnetic 
moment is subtracted from all vertices in these diagrams. Technically the gauge invariant 
anomalous magnetic moment is the hardest entry in the expression for the vertex at small 
transferred momenta, and this prompts its separation. It turns out that the anomalous 
magnetic moment does not generate radiative-recoil corrections (see, e.g., [21 II])- 




FIG. 1: Electron polarization insertions 




FIG. 2: Muon polarization insertions 

The scale of integration momenta in the graphs in Fig. [T] and Fig. [2] is determined by the 
muon mass, and is much higher than external virtualities determined by the characteristic 
atomic momenta. In order to obtain radiative-recoil corrections it is sufficient to calculate 
matrix elements of the diagrams in Fig. [T] and Fig. [2] with the mass-shell external momenta 
with vanishing spatial components. Then contribution to hyperfme splitting is given by the 
sum of matrix elements of the gauge invariant sets of diagrams in Fig. [T] and Fig. [2] calculated 
between the free electron and muon spinors and multiplied by the square of the Coulomb- 
Schrodinger wave function at the origin. Explicitly contributions to HFS generated by each 
of the gauge invariant sets of graphs in Fig. [T] and Fig. [2] have the form 

( 1 1 \ Za m 

X \k 2 + 4fik + i0 + kz-4fik + i0 ) = ~k~ ~~M ' 

where m and M are the electron and muon masses, respectively, fi = m/(2M), k is dimen- 
sionless exchange photon momentum measured in terms of the muon mass, a is the fine 
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structure constant, Ep = (8/3)(Za) 4 (m/M)m is the Fermi hyperfine splitting energy, L^ u 
is the muon line factor, and the broken brackets denote matrix elements between hyperfine 
states. Constant Z measures the muon charge in terms of the electron charge and is equal 
one in muonium, but we keep it as an indicator on the origin of different contributions. 




FIG. 3: Two-photon exchange 



Generically dimensionless recoil contributions Ae are enhanced by powers of large log- 
arithm of the muon-electron mass ratio (see, e.g, [HE])- Origin of these large logarithms 
is quite transparent. In the case of the skeleton muon factor momentum integral in 
Eq. is linearly divergent at low integration momenta of order mZa (we temporarily re- 
turn to dimensionful momenta in these considerations). This linear divergence corresponds 
to the classic nonrecoil Fermi contribution to HFS and should be subtracted in calculation 
of recoil contribution. After subtraction the skeleton integral in Eq. ^ becomes logarithmic 
in the wide integration region m < k < M and generates the leading recoil correction to 
HFS. Considering one-loop radiative insertion in the muon line in Fig. |4]we need to separate 
gauge invariant one-loop anomalous magnetic moment and the remaining muon factor. The 
one-loop anomalous magnetic moment has the same momentum dependence as the skeleton 
contribution, and generates only nonrecoil contribution. Due to generalized low-energy the- 
orem for virtual Compton scattering (31 [5] the remaining one-loop muon factor is suppressed 
by an additional factor k 2 /M 2 in comparison with the skeleton muon factor. This suppres- 
sion not only makes the integral in Eq. convergent in the infrared region, but also makes 
the integral nonlogarithmic. As a result one-loop muon factor generates only nonlogarithmic 
radiative-recoil corrections obtained in [21 IH [5, 6J. Insertion of vacuum polarization in the 
radiative photons does not change the low momenta behavior of the muon factor and 
the integral over exchanged momenta in Eq. corresponding to the diagrams in Fig. [I] and 
Fig. [2] remains convergent and nonlogarithmic. Then the leading recoil correction to HFS 
generated by the diagrams in Fig. [2] with the muon polarization insertions does not contain 
large logarithms and is a pure number. The case of of electron polarization insertions in 
Fig. [l] is a bit more complicated. While it is still true that integration over the exchanged 
momentum is nonlogarithmic, another source of large logarithms arises in these diagrams. 
As we already discussed all characteristic momenta in these diagrams are or order of the 
muon mass M . Then the electron polarization insertions in the diagrams in Fig. [T] enter 
in the asymptotic regime, and the leading contribution to HFS generated by the diagrams 
with the electron polarization insertions in Fig. [l] is just the product of the leading asymp- 
totic term in the high momentum expansion of the electron polarization operator and the 
radiative-recoil correction to HFS generated by the one-loop muon factor. In other words 
we can say that the leading logarithmic contribution to HFS generated by the diagrams in 
Fig. [T] is obtained from the respective nonlogarithmic contribution of the diagrams in Fig. [4] 
by substitution of the running coupling constant a(M) for radiative photons. 
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FIG. 4: One-loop radiative insertions in the muon line 

II. ANALYTIC CALCULATIONS 

Actual calculations start with consideration of separate contributions generated by the 
diagrams with self-energy, vertex, and spanning photon insertions in Fig. [T] and Fig. [2j 
Redefining the muon factor by extracting from it all factors of a, we write the general 
expression in Eq. ([I]) in the form 

AE = Ae a{Z2 f Za) ^E F , (2) 

where 

Ae = Ae s + 2Ae A + Ae s . (3) 

Explicit expressions for different terms on the right hand side in Eq. ^ are derived in the 
Feynman gauge. We start with the respective diagrams in Fig. [4] without polarization inser- 
tions but with massive radiative (but not exchanged) photons. We account for polarization 
insertions by considering photon mass squared to be A 2 = 16yU 2 / (1 — v 2 ) or A 2 = 4/(1 — v 2 ), 
for the electron and muon polarizations, respectively. Then insertion of the polarization 
operator means an additional integration over velocity v with the weight 

dv 9 / v 2 \ 

v 2 l-^r • (4) 



1 - v 2 V 3 , 

Below we collect integral representations for the contributions Ae^ obtained directly from 
the Feynman diagrams. For the contribution of the self-energy diagram we obtained 



Ax 



2vr 2 J J J 1 - v 2 V 3 J J k A k A - 16fx 2 k 2 

Ae S i + Ae E2 , 



x<i3^o L r^-(3^-2fc 2 )(l-x) 2 ^ 2(1_K '' )/y 



A r 



1 

aT 



where 



A (x) = x 2 + A 2 (l - x) = y(l - x)a\(x, y), 
Ax(x,y) = y(l - x){-k 2 + 2k + a\ - iO). 
The spanning photon contribution has the form 



i f 1 , f X , < , f 1 dv 2 / v 2 \ fd 4 k 2k 2 
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l-2y -2 + 2x(l - y) 



-A 2 (l-x)y 



A 



—6bk 
where 
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2y 2x(l -y)+x 2 



X 2 (l 



A 2 
x){2-y) 



A 2 



A 2 A 2 

A(x,y) = y(l-y)(-k 2 
x 2 + A 2 (l - x) 



+ 



2bk y{\ - y) 
A 2 



= Ae=i + Ae-. 



a (x,y) 



1/(1-1/) 



2bk + a 2 -iO), 
&(x,y) = 1 



The vertex contribution is as follows 

• i 
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2tt 2 
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Ae Ai + Ae A 2 + Ae A3 , 



where 



A u (x,y,u) = y(l - y){-k 2 + 2b u k + a 2 ), 

2 , N v r __ s 1 — y — xu + yu 



a u( x , V, u ) = ua (x, y), K(x, y, u) 



i-y 



Our next goal is to calculate the momenta integrals in Eq. Eq. ([6]), and Eq. (J7|) 
analytically. We explained in Section [I] that due to soft low-energy behavior of the muon 
factor L^ v the momentum integral in Eq. ([I]) is convergent and nonlogarithmic. Hence, 
Ae does not depend on the electron-muon ratio explicitly (it depends only on the effective 
photon mass A), and one can safely omit factor [i in the electron denominators in the 
square brackets in Eq. 0, what significantly simplifies all calculations. However, each of 
the exchanged momentum integrals in Eq. ([5]), Eq. (|6]), and Eq. ^ corresponding to the 
self-energy, the spanning photon, and the vertex insertion contributions to the muon factor 
is logarithmic in the wide integration momenta region m/M < k < 1. At small integration 
momenta of order m/M these would be divergences are cut off by the /^-dependent electron 
denominators. This is the reason why we preserved parameter fi in sums of the electron 
propagators 2k 2 / (k 4 — lQ^k^) in Eq. (|5]), Eq. ([6]), and Eq. Q. Separate contributions in 
Eq. ([5]), Eq. (|6]), and Eq. §f§ are gauge noninvariant, and one could hope to find a gauge 
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where individual contributions are infrared finite even if one neglects the electron mass in 
the electron denominators. This is not the case, such gauge does not exist and infrared 
divergences survive in individual contributions even in the infrared soft Yennie gauge [3]. 

We would like to get rid of the would be logarithmic divergencies at the level of integrands, 
so that divergent for \x = momentum integrals would not arise at all. To this end we 
wrote all integrals in Eq. (|5]), Eq. ([6]), and Eq. ^ in terms of similar Feynman parameters 
and organized them in three groups, combining contributions from different integrals. We 
obtained the representation for the energy splitting in the form 



i=3 



Ae = Ae s + 2 Ae A + Ae s = ^ Ae i 



(8) 



i=l 



where 



Aei = Ae S2 + Ae 31 + 2Ae A i, (9) 

Ae 2 = Ae E i + Ae H2 + 2Ae A2 , (10) 

Ae 3 = 2A eA3 . (11) 

Momentum integrals in Eq. (J9| and Eq. (10) are apparently logarithmically divergent at 
[i = 0, and require additional transformations. 

Consider first Aei in Eq. ([9]). After some integrations by parts over the Feynman param- 
eters it acquires the form 



A ei = 
where 



7T- 



1 fX rl 

dx dy 
!i Jo Jo 



dv 



l-v< 



VI i 



v 
~3 



d 4 k 



k 2 (k 4 - 16 fi 2 k 2 ) 



2 Ay 



2J A 1 k 2 B 1 
A A 2 



(12) 



..-I j (.r. //. A) = 3( l- 2 ^- ) -.V- :ir 



, Ax \ o A 2 (l — x)v , 
lr)u+ { 2-— jy 2 + jy (-9 + 3x-2y), 



and 



B l 



xl 1 



X 



2 x^ x(l-x)y 



Ar 



y(l-y) + (-2 + x) 
y 2 (x - y) 



xy 
A^ 



+ [-2 + 2x(l-y)-x 2 -X 2 (l-x)y] 

After the Wick rotation the integral for Aei i n four-dimensional Euclidean spherical coor- 
dinates (k = kcos6, \k\ — ksinO) has the form 



Aei 



dx / dy 



1 dv 



1-v 2 



:V 2 1 



k 2 dk 2 



X 



At k 2 B x d 

+ 



3 J J k 4 + Wn 2 k 2 cos6 2 n J() 
k 2 + a 2 



P7V 

/ d9 sin 2 6(2 + cos 2 
Jo 



(13) 



2/(1-3/) y 2 (l-y) 2 da 2 



(k 2 + a 2 ) 2 + Ab 2 k 2 cos 2 9 
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Now we see that omitting the term with /i 2 in the denominator would result in an apparently 
logarithmically divergent integral. However, this is not the case since the coefficient before 
the would be divergence turns into zero due to the identity 



dx 



dy 



An 



0. 



(14) 



that is valid for any A. This means that the momentum integral in Eq. (13) is convergent 
even at \i = 0. We obtain an explicitly convergent form for this integral by subtraction from 
the integrand taken at /x = its leading low momentum asymptotic term. This subtraction 
does not change the integral since the coefficient before the the subtracted term is identically 
zero due to the identity in Eq. (14). Of course, we have chosen this particular combination of 
terms on the right hand side in Eq. ^ precisely because we wanted to organize an infrared 
finite contribution. After subtraction and momentum integration we obtain a finite integral 
representation for Aei 



Ad 



dx 



3 , 



dy 

3 1-y 
Ao + y(l- x y 



dv 



1-v 2 
(2L + Li) 



v 
~3 



9 A x 1 

: — hi — ^ 

4 An a 2 



2y 2 (l-x) 2 



(L + L x ) 



where L n = b 2 dss n / (a 2 + b 2 s) } or explicitly 

a 2 + b 2 



A) 



In 



L 1 = l 



a 
V 2 



Lq. 



Consider now the contribution in Eq. (10) that is defined by the integral 

•i 



Ae 2 
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Xfcn 
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dx 
x(l — x 



dy 



dv 
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1 -t» 2 



1 - 



d A k 



2k 2 



_ A k 2 x(l - x)y(l - y) 
A AA A 2 



k 4 k 4 - 16/x 2 fc 2 

2 



k 2 B, 



where 



(15) 



(16) 



(17) 



A 2 {x,y, A) 

,2 



[1 - x)(2x - y) + (-2 + x 



2 

XV 



B 2 = y 3x - 2y 



x 
~2 



3x 2 y . 2 
~~ ^ — I" x u 



+ (-2 + x) 1 + 



3x 



y 



x 2 y 2 



After the Wick rotation we obtain in four-dimensional Euclidean spherical coordinates an 
integral 



x 



Ae 2 



d6 sin 2 9 cos 2 9 



dx 



dy 
A 2 



dv 



+ 



x(l 



1 



V 



2b 



k 4 dk 2 



k 4 + lQp?k 2 cos9 2 
d 1 1 



(18) 



Ti j i y{l-y)k 2 ' A y 2 {l - y) 2 da 2 J (k 2 + a 2 ) 2 + 4b 2 k 2 cos 2 0' 

that is apparently logarithmically divergent for \x = 0. Like in the case of the contribution 
Aei the coefficient before the would be divergence turns into zero due to the identity 
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1 r (l-x)y 
dx \ dy — ^ — A 2 (x,y, A) = 0, 



o 



A 



(19) 



o 



that is valid for any A. Using this identity we subtract the would be divergent at \i = 
integral of the leading low momentum asymptotic term in the integrand in Eq. (18), and 
obtain a finite integral representation for Ae 2 



Ae 5 



dx 
A 2 



2bA 



dy 



dv 



1-v 2 
x(l - y) 



v 2 \ 1 



V 

~3 



[1 - x)yA 2 



<2A\ 



In — + L 

a 2 



(20) 
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(L - Li) + 



Bo 



y{l-x)A 



Li 



The momentum integral in Eq. (11) is convergent even at jj, — 0, and since we are 



interested in the leading recoil contribution we calculate it in this limit 



xj (L 0u + L lu ) 
where 



Ae, 



(2 - 3x) 



dx 

y(x - y) 
A 



dy 



dv 



vO- -y) Jo l ~v 2 

-[l-x + (x- y) 2 ] 



v 2 1- 
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o K 



+ (L 0n - 3L lu )x(l - X 
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An 



T - 1 U T 

Of, 



(22) 



III. NUMERICAL CALCULATIONS AND RESULTS 



Let us notice that the expressions for Ae, obtained in Eq. (15), Eq. (20), and Eq. (21) 



allow us to reproduce the radiative-recoil correction to HFS generated by the diagrams in 
Fig. [4] with one-loop radiative insertions in the muon line. To obtain these corrections we 
need to omit integration over velocity v and all explicit factors containing v in Eq. (15), 
Eq. (20), and Eq. (21), and also let the radiative photon mass A to be zero. The resulting 



integrals admit analytic calculation and the results are collected in Table [T| Restoring the 
overall dimensional factor we obtain 



AE = Ae^ q -—Ep = -C 3 - 3vr 2 In 2 + 

7T Z M \Z 



39\ {Z 2 a)(Za) m 



7P 



M 



(23) 



what precisely reproduces the old result [31 E] ■ 

Now we are ready to proceed to calculations of the contributions to HFS generated by the 
diagrams in Fig. [T] and Fig. [2] The leading logarithmic term corresponding to the diagrams 
in Fig. [T] with the electron polarization insertions in the radiative photon is obtained by 
substituting the running coupling constant a(M) for the radiative photon in the result in 
Eq. (23), and is equal to (we do not include a in the definition of Ae, it will be restored in 
the finite results) 
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TABLE I: One-loop radiative recoil corrections 



Ae[ 0) -15C(3)+6^ 2 ln2-5^ 2 + 5f 
A4 0) fC(3)- 3f]n2 + £-! 

Aef l|i C (3)_ 3^i n2 + l^-24 



Ac(°) |C(3)-37r 2 ln2 + f 



A^= 2 ln-A e (°>. 
3 m 



(24) 



Similar equations are also valid for each of the contributions Ae- ''. For practical calculations 
we defined nonlogarithmic contributions to HFS as 



AC, = Ae 



3 m 9 



(25) 



We calculated contributions in Eq. (25) numerically using the integral representations in 
Eq. (15), Eq. (20), and Eq. (21), with the photon mass squared A 2 = 16/z 2 /(l — v 2 ), and 



obtained 



Ad = -1.4965(7), AC 2 = 2.2442(6), AC 3 = 5.7905(5). 



(26) 



Then the total contribution to HFS due to electron polarization insertions in Fig. [I] has the 
form 



Ae e = ( 3C(3) - 2tt 2 In 2 + — J In — + 12.227(2). 



4 J 



m 



(27) 



The contribution to HFS due to muon polarization insertions in Fig. [2] was calculated 



numerically with the help of the integral representations in Eq. (15), Eq. ^201) , and Eq. (21). 



This time A 2 = 4/(1 — v 2 ) plays the role of the photon mass squared in the integrands, and 
we obtained 



Aei = 0.2987, Ae 2 = -0.0003, Ae 3 



-1.2291, 



or 



(28) 



Ae, 



-0.9307. 



Restoring all dimensional factors we obtain from Eq. (27) and Eq. (29) 

a(Z 2 a)(Za) m 



AE P 



3C(3) - 2vr 2 In 2 + — ) In — + 12.227(2) 
4 / m 



71'' 



M~ 



AE„ 



-0.931 ^ Za) ^E F . 



(29) 



(30) 



(31) 



7T 3 M 

Total contribution to HFS due to diagrams with electron and muon polarization insertions 
in Fig. [T] and Fig. [2] is 



9 



AE 



3C(3) - 27r 2 ln2 + 



13 



M 

In — + 11.297(2) 

m 



a 3 m 
^1 



E 



F- 



(32) 



This result together with the results [TJ, [HI El HO] completes calculation of all radiative- 
recoil corrections of order a 3 (m/M)Ep generated by the diagrams with electron and muon 
polarization insertions. 
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